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The electronic and magnetic properties of the mother material LaOFeAs of new super- 
conductors have been carefully studied using first-principles electronic structure calculations 



Oh: 

i^ . based on the generalized gradient approximation in the density functional theory. The present 

j^ . calculation predicts that the ground state of LaOFeAs is antiferromagnetic with a stripe type 
magnetic moment alignment leading to orthorhombic symmetry of the crystal. In this par- 

1^ . ticular magnetic state, the density of states at the Fermi level is very small. On the other 

S ■ hand, LaOFeP has turned out to be paramagnetic and a good metal. Implications of the 

O . results regarding the experimental observations will also be presented. 
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Discovery of a new superconductor always attracts strong attention of the scientific com- 

• ^H , 

^ ■ munity particularly if the critical temperature Tc is above 20 K. It is rather surprising that 

C^ ' immediately after the recent news of superconductivity of La0i_2,.F2.FeAs^ several papers, 

both theoretical and experimental, have been circulated in the community. There may be four 
main reasons why this discovery has produced such a strong impact. First, Tc is relatively 
high 26 K. Second, it is a new finding that Fe, which is a typical magnetic element, seems 
to be participating to superconductivity. Third, the crystal takes again a layered structure 
and the doping region and the superconducting region are geometrically and electronically 
separated. Fourth, the material seems to have strong flexibility in the choice of constituent 
elements suggesting possibility of higher Tc materials in the same category.^ 

In order to go further in exploring possibility of higher Tc materials, it is essential to 
understand the basic electronic properties of the mother material LaOFeAs and the role of 
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partial replacement of O with F. In the analogy of high T^ cuprates,^ important questions 
about LaOFeAs may be whether magnetism is involved or not and whether it is metallic or 
insulating. These questions have already been addressed by other works^"^ in which standard 
DFT band calculations have been performed. At least within our knowledge, it was concluded 
that the system is metallic and nonmagnetic with possible strong antiferromagnetic (AFM) 
fluctuation. We have also performed similar calculations and have arrived at a conclusion that 
the ground state will be a particular AFM state. The density of states (DOS) at the Fermi 
level is very low suggesting that the system may be a bad metal. 

Briefly, our calculation is based on the PAW method^ and the PBE version^ of generalized 
gradient approximation in the density functional theory. Spin-orbit interaction was not taken 
into account. We used our in-house computer code named QMAS (Quantum MAterials Sim- 
ulator) and have tested the reliability of the code in several ways. In relation to the present 
system, we have performed test calculations for LaOMnAs, LaONiP and MnAs with NiAs 
structure. In all of these test calculations, we have confirmed consistency between our results 
and other existing ones^'^ or experimental facts. ^'^ In the actual calculations for LaOFeAs, 
we have studied six different magnetic states: 1) nonmagnetic (NM), 2) ferromagnetic (FM), 
3) AFM with G-type (AFM-G, Fig. la), 4) AFM with S-type (AFM-S, Fig. lb), 5) AFM with 
CEl Type (AFM-CEl, Fig.lc) and 6) AFM with CE2 type (AFM-CE2, Fig.ld). 

The convergence criteria were 10~^ electrons/bohr'^ for the self-consistent-charge distance, 
5 X 10^^ hartree/bohr for the maximum force amplitude and 5 x 10^^ hartree/bohr^ for the sum 
of diagonal components of the stress tensor, respectively. The plane- wave energy cutoff was set 
to 20 hartree. The number of k points was set to 12x12x6 for the 1x1 lattice (NM, FM and 
AFM-G), 10x10x6 for the c(2x2) lattice (AFM-S), 6x12x6 for the 2x1 lattice (AFM-CEl), 
or 6x6x6 for the 2x2 lattice (AFM-CE2). 

In order to make a guess which magnetic state is the most stable among the above- 
mentioned six candidates, we first estimated the total energy of each magnetic state using 
the common crystal structure determined experimentally.^^ The results are shown in the 2nd 
column of Table I as the difference from the total energy for the corresponding NM state 
/S.E = Ef^y^-jpj^^ — Enmi where E/^^-jp]^ and EjyM are total energies for (A)FM and NM states. 
The magnetic moment of Fe estimated within the radius of 2.2 bohr is shown in the 3rd 
column. Thus, the AFM-S state is the most stable among the above six candidates (NM, FM, 
AFM-G, AFM-S, AFM-CEl and AFM-CE2). Knowing the most stable magnetic state at this 
level, we optimized the crystal structure (lattice constants and internal coordinates) for NM 
and AFM-S states to evaluate the AFM stability energy more precisely. We also optimized 
the crystal structure for AFM-G because of the reasons mentioned below. The revised total 
energy differences and magnetic moments after structural optimization are shown in the 4th 
and 5th columns, respectively, of Table I. The magnetic stability energy for both AFM-G 
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and AFM-S is reduced because the effect of structural optimization is more significant for the 
NM state. In the following, we will discuss basic features of the electronic structures of some 
selected magnetic states. 

Figure 2 shows the density of states (DOS) for the NM state of LaOFeAs with the ex- 
perimental structure/^ which is very similar to the corresponding one in other works. ^'^ The 
whole valence states are divided into p bands coming from O and As and d bands coming 
from Fe. 

The width of d band is about 4 eV, being about 2/3 of that of metallic bcc Fe and 
suggesting significant itinerant character of d states. By counting the capacity of bands and 
number of valence electrons, the formal number of d electrons in LaOFeAs is six per formula 
unit and the Fermi level is located at the shoulder just below the dip in the d band. We have 
found that there is only very little contribution from oxygen p states to the states at the Fermi 
level. As is well known in the physics of transition metals, ^^ this situation of the Fermi level 
with regard to the d band may favor AFM instability rather than FM one if any magnetic 
instability may exist. This is indeed the case as Table I shows. However, it is fairly difficult 
to find the most stable AFM states because there are infinite number of possible metastable 
AFM states. In the present system, as Fe forms a square lattice, the AFM-G state may naively 
be expected to be realized. With the experimental crystal structure, the AFM-G state is more 
stable than the NM state by 0.08 eV per formula unit. (Note that the AFM-G stability energy 
with respect to the NM state is reduced to 0.016 eV after structural optimizatin. See Tablel.) 
However, this AFM state has a rather strange feature in DOS, which is shown in Fig. 3. Due to 
the G-type arrangement of the exchange potential, d states will have a tendency of localization. 
This tendency is particularly strong for d322_^2 due to both the 2 dimensional network and the 
flattening of the tetrahedral coordination of As. This strong localization of d-^2'2_j.2 produces 
a rather sharp peak in the region of a dip in DOS where the Fermi level is located. With 
such a sharp peak in DOS at the Fermi level, we generally expect some structural distortion. 
Although we have tested several possible structural distortion models, none of them could 
remove the sharp peak from the Fermi level. 

Instead of exploring further possibilities of structural distortion within the AFM-G mag- 
netic structure, we studied three other types of AFM spin arrangement, which reduce the 
symmetry of the crystal from tetragonal (for NM, FM and AFM-G) to orthorhombic. These 
magnetic structures look strange at first sight, but the angular variation of d orbitals can 
adjust the covalency between Fe atoms to be consistent with the magnetic structure. In other 
words, magnetic structure is accompanied with orbital ordering and then crystal structure is 
consistently modified from the tetragonal symmetry. Such examples are frequently observed 
in perovskite transition metal oxides. ^^ Our calculation predicts that the AFM-S state is 
the most stable among six possible magnetic states studied. We therefore present some more 
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details of the optimized crystal structure and the electronic structure of the AFM-S state 
below. 

As mentioned above, the orthorhombic symmetry of the magnetic moment arrangement 
in the AFM-S state will lead to modification of the crystal structure mediated by orbital 
ordering, which can be seen in the spin density distribution in Fig. 4. Qualitatively, the states 
of local majority spin state with energy near the Fermi level contribute to the spin density. 
As the states in this energy range are more or less antibonding between nearest neighbor Fe 
atoms, the spin density is more localized along the stripe of FM spin alignment. According to 
the structural optimization, the lattice constants of the orthorhombic cell, which is originally 
the c(2x2) cell (a = 4.03533 x ^/2 = 5.70682 A, c = 8.74090 A) of the tetragonal lattice, are 
a = 5.686 A (parallel to the stripe), b = 5.758 A (perpendicular to the stripe) and c = 8.707 
A. The z values for La and As {zLa and zas) are 0.142 and 0.646, respectively. For comparison, 
the optimized tetragonal lattice constants for the NM state are a = 5.677 A, c = 8.640 A 
with ZLa = 0.145 and zas = 0.639. Figure 5 shows the total DOS and Fe d partial DOS of the 
AFM-S state. Orbital decomposed partial DOSs in the AFM-S state are shown in Fig. 6 The 
band dispersion in the Brillouin zone of the c(2x2) unit cell of the original lattice is shown in 
Fig. 7. Although a couple of bands cross the Fermi level, very small DOS at the Fermi level 
(0.246 states/eV/cell for each spin) implies that the system may be a bad metal. The states 
at the Fermi level have contributions dominantly from minority-spin dyz and d3^2_^2 orbitals 
and less dominantly from majority-spin dyz and d^2_j^2 orbitals. The contribution from As 
Py is noticeable but an order of magnitude smaller than the contribution from Fe d orbitals. 
Note that the x (y) axis is parallel (perpendicular) to the stripe of the AFM-S structure. 

As for the electronic structure of the NM LaOFeP, we obtained very similar results to 
those reported by Kamihara et alM Now we ask ourselves how different it is between LaOFeAs 
and LaOFeP in our band structure calculations. Their NM states have nearly identical DOS 
except that the band width is slightly narrower in LaOFeAs than in LaOFeP. There is also 
some noticeable difference in the contributions from As and P to the states at the Fermi level: 
p state contribution from As is much larger than that from P. Such a subtle difference in the 
NM state produces striking difference between them in the magnetic state. We have found 
that the AFM-S state converges to the NM state very slowly in the process of structural 
optimization for LaOFeP. Therefore, LaOFeP is paramagnetic and a good metal with possible 
AFM fluctuation corresponding to the S type. 

The results of the present calculation are qualitatively consistent with some experimen- 
tal data such as resistivity.^' ^^ As was already mentioned above, the calculated electronic 
structure implies that LaOFeAs may be a bad metal while LaOFeP may be a good metal. 
Qualitatively such behavior is actually observed experimentally.^' ^^ Stability of AFM state 
and very small DOS at the Fermi level for LaOFeAs may be consistent with the fact that the 



4/6 



J. Phys. Soc. Jpn. Letter 

undoped pure system does not show superconductivity. Note that pure LaOFeP is a supercon- 
ductor below 4-7 K.^^ On the other hand, we do not know at the present stage how to explain 
the observed magnetic susceptibility of LaOFeAs^ in terms of AFM ordering and small DOS 
at the Fermi level. Experimentally, the susceptibility is almost temperature independent like 
a metallic Pauli susceptibility and the magnitude of the order of about 0.5 x 10~^ emu/mol is 
too large to be explained by DOS at the Fermi level. The role of partial oxygen replacement 
with F is now being studied. 

We thank Prof. H. Fukuyama for useful discussion and Dr. M. Kohyama, Dr. S. Tanaka 
and Dr. T. Tamura for their help in developing our computational code QMAS. The present 
work is partially supported by the Next Generation Supercomputer Project, Nanoscience 
Program and also partly by Grant-in- Aids for Scientific Research in Priority Area "Anomalous 
Quantum Materials", both from MEXT, Japan. The calculations were performed using the 
AIST Super Cluster at the Tsukuba Advanced Computing Center (TACC), AIST. 
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Fig. 1. Antiferromagnetic order patterns investigated in this study, a: AFM-G, b: AFM-S, c: AFM- 
CEl, d: AFM-CE2. Squares and a rectangular in broken lines represent unit cells for the corre- 
sponding AFM states. Schematic view of the experimetal unit cell is shown also. 



Fig. 2. Density of states (DOS) for the NM state of LaOFeAs with the experimental structure. The 
energy zero corresponds to the Fermi level. 



Fig. 3. Density of states (DOS) for the AFM-G state of LaOFeAs with the experimental structure. 
The energy zero corresponds to the Fermi level. The upper half represents the DOS for up spin 
while the lower half is for down spin. 
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Fig. 4. Spin density distribution viewed along the c direction for the AFM-S state of LaOFeAs. 
Isodensity surface contours for the difference between the up-spin and down-spin electron densities 
are plotted. The blue and red surfaces correspond to ± 0.005 electrons/bohr"^, respectively. Larger 
surfaces are located on Fe atoms. As atoms arc at {x,y) = (0.25,0.25), (0.75,0.25), (0.25,0.75) and 
(0.75,0.75). The boundary corresponds to the square shown in Fig. lb. Along a (b), the lattice 
shrinks (expands) from the original square. 



Fig. 5. Total DOS (a) and Fc d partial DOS (b) for the AFM-S state of LaOFeAs. The energy zero 
corresponds to the Fermi level. The upper half represents the DOS for up spin while the lower 
half is for down spin. 



Fig. 6. Orbital decomposed partial DOS for the AFM-S state of LaOFcAs. The energy zero corre- 
sponds to the Fermi level. The upper half represents the DOS for up spin while the lower half is 
for down spin. 



Fig. 7. Band dispersion in the Brillouin zone of the c(2 x 2) unit cell of the original lattice for the AFM- 
S state of LaOFeAs. The right panel is a magnified drawing of the left one near the Fermi level. 
The energy zero corresponds to the Fermi level. Wave vector is shown in the unit of (7r/a,7r/6,7r/c). 



Table L 2nd column: Total energy of each magnetic state with reference to that of NM state with the 
common lattice structure from experiment, ^^ 3rd column: Magnetic moment of Fe, 4th column: 
the same quantity as that of column 2 after structural optimization only for AFM-G and AFM-S, 
5th column: Magnetic moment of Fe after structural optimization. 

Type A^* (eV/formula) M** (^b) AEopt* (eV/formula) M*;^ (^b) 

FM 

AFM-G -0.080 1.89 -0.016 1.51 

AFM-S -0.174 2.12 -0.098 2.00 

AFM-CEl 

AFM-CE2 

Enm 
Magnetic moments were estimated within a radius of 2.2 bohr. 
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